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For collectively interacting repulsive particles driven on triangular substrates, we show that for 
certain directions of drive a spontaneous symmetry breaking phenomena occurs where the particles 
can flow in one of two directions that are not aligned with the external drive, giving rise to a positive 
or negative Hall current. Along these directions, the particle flow is highly ordered, while in the 
direction of the drive the flow is disordered. We also find a number of dynamical phase transitions 
and unusual hysteretic properties that arise due to the symmetry breaking properties of the flows. 
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A number of systems can be represented as collectively 
interacting particles moving over an underlying periodic 
substrate, such as models of atomic friction [1], vortex 
flow in superconductors with periodic arrays of pinning 
sites [2-4], and colloids interacting with periodic opti- 
cal trap arrays [5-7]. Typically these systems exhibit a 
rich variety of dynamical phase transitions under increas- 
ing applied drive, such as transitions from highly disor- 
dered fluctuating flows to highly ordered one-dimensional 
(ID) and 2D flows, with strong hysteresis between the 
phases [1,3,4]. Other effects such as negative resistance 
and switching can also occur in conjunction with these 
dynamic transitions [1,3]. Usually in these systems the 
average drift velocity is finite in the direction of the ap- 
plied force but zero in the direction transverse to the 
force. Under certain conditions, however, the particles 
do not flow in the direction of drive, but lock into a sym- 
metry direction of the underlying lattice. For a square 
substrate lattice, changing the angle of drive with respect 
to the substrate produces a series of dynamical locking 
phases forming a Devil's staircase structure [4]. Lock- 
ing effects have been demonstrated for vortices moving 
in superconductors with periodic pinning arrays in both 
simulations [4] and experiments [8] , theoretically for elec- 
trons on periodic substrates [9], and in experiments [5] 
and theory [10] for colloids driven over periodic traps, 
where it was shown that the locking can be used to sort 
particles [11]. The directional locking effects are not due 
to collective interactions of the driven particles and can 
occur for a single driven particle. 

For certain periodic substrate lattice geometries, there 
are two symmetry directions of the substrate that the 
particles could follow instead of the direction of drive. 
An example is a triangular lattice oriented such that the 
applied drive is between two symmetry directions, as il- 
lustrated schematically in Fig. 1. In this work we show 
that when there are many interacting driven particles in 
such a geometry, spontaneous symmetry breaking can oc- 
cur where the particles preferentially flow globally along 
either of the two symmetry directions rather than in the 
direction of drive. Since the substrate is symmetric, the 
occurrence of the global flow can be viewed as a spon- 



taneous symmetry breaking. This transition is driven 
by collective interactions between the repulsive particles. 
In the limit of a single particle, symmetry breaking is 
not observed. When the global flow follows one of the 
symmetry directions of the substrate, the flow is ordered 
with few fluctuations and a uniform spacing between the 
particles, reducing their interaction energy. If the global 
flow instead follows the direction of the drive, the par- 
ticle configuration and fiow are much more disordered. 
At T = for ordered configurations of the driven parti- 
cles, such as at commensurate fillings, the particles fiow 
in the direction of drive. At incommensurate fillings or 
finite temperatures, however, enough perturbations can 
occur to allow the system to order into one of the two easy 
fiow directions. For different system parameters such as 
waiting times, ramp rates, fillings, and initial conditions, 




FIG. 1. Schematic of the system with a triangular sub- 
strate. The potential minima are at the center of the circles. 
Each minima captures one or more interacting particles. The 
applied driving force is fdy (solid arrow). There are two easy 
flow directions at ±30° from the drive (dashed arrows). 
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the symmetry breaking is equally likely to occur in either 
direction; thus, for an ensemble of realizations, the sym- 
metry is restored in the statistical sense. Additionally we 
find a series of dynamical phase transitions where the par- 
ticle flow changes from ordered to disordered. This can 
produce unusual hysteretic effects. For example, when in- 
creasing the drive, the transverse velocity may lock to the 
positive direction at low drives, but if the system passes 
through a disordered flow regime at high drives, then 
when decreasing the drive the flow has an equal prob- 
ability to lock in the negative direction. We explicitly 
show this effect for colloids moving over periodic sub- 
strates and vortices in superconductors with triangular 
pinning arrays. 

We consider a system of Nc overdamped repulsively 
interacting particles interacting with an underlying tri- 
angular periodic substrate with N minima. The equation 
of motion for a particle i is 

^ = f,, + f, + f, + . (1) 

Here the particle-particle interaction force is fy = 

~ ^^^(^ijOj where we use a Yukawa or screened 
Coulomb interaction potential appropriate for colloids, 
given by V{rij) = (<?V|ri - r^l) exp(-K|ri - rj|). Here q 
is the colloid charge, 1 /k is the screening length, and r^jj) 
is the position of particle We also consider vortices 
in thin film superconductors with periodic pinning. In 
this case the vortices interact with a logarithmic potential 
V{rij) = — In(rij). For the colloids we consider a trian- 
gular substrate which can capture multiple colloids per 
minima: fs = X^^'^^ Asin(27rpi/ao)[cos(6'i)x - sin(6'i)y], 
where pi = xcos{6i) — ysm{0i) + ao/2, 9i = 7r/6, 
O2 = 7r/2, and 63 = 5tt/6. Such substrates can be cre- 
ated with optical trap arrays [7] and have been studied in 
previous simulations [12]. In the case of vortices, the pin- 
ning sites are modeled as short range attractive parabolic 
pinning sites in a triangular lattice that capture one vor- 
tex each. The additional vortices sit in the interstitial 
regions between the pinned vortices. Length is measured 
in units of the substrate lattice constant ao and we use a 
screening length for the colloids of k = 3/ao. The ther- 
mal force fr is a randomly fluctuating force from random 
kicks. We consider both T = and finite T. The initial 
particle positions are obtainc;d by simulated annealing 
from a high temperature. The driving force fd is applied 
in the y direction as shown in Fig. 1. The drive is in- 
creased in small increments to a finite value and then 
similarly decreased back to zero. We measure both the 
velocity in the direction of drive, Vy = J^i' -y, and the 
velocity transverse to the drive, 14 = Vy • x. 

We first consider the case of colloidal particles inter- 
acting with a triangular substrate. We concentrate on a 
regime where the number of colloids is greater than the 
number of potential minima, as in recent experiments on 
triangular substrates [7] . Fig. 1 shows a schematic of 




FIG. 2. (a) Velocity in the direction of drive Vy vs applied 
drive fa for a system of interacting colloids on a triangular 
substrate for Nc/N = 1.39. (b) Corresponding transverse 
velocity Vx ys fd. Regions I-IV correspond to different particle 
flow regimes, (c) The transverse velocity Vx vs fd for the same 
system as above but for slightly different initial conditions 
showing that the transverse velocity reverses from (b) . 

our system. The dashed arrows indicate the symmetry 
directions for particle flow. In Fig. 2(a) we plot the av- 
erage particle velocity in the direction of the drive, Vy, 
vs applied force for a system with Nc/N = 1.39. In 
Fig. 2(b) we show the corresponding average velocity in 
the x-direction, 14 . Here four distinct dynamic regions 
can be observed. The flrst region, I, is the pinned phase 
where particles do not move in either direction. Region 
II is the spontaneous symmetry breaking phase where 
the global flow of particles follows either the positive or 
negative 30° direction as in Fig. 1. In this regime Vy and 
Vx both increase linearly with increasing fd- In Fig. 2(b) 
the particles flow along —30°, producing a negative "I4 
or Hall current. As fd increases, the particles start to 
move in the direction of the drive and there is a transi- 
tion to region III flow. This transition appears as a sharp 
drop in I4 with a corresponding increase in the slope of 
Vy. Region HI is a disordered flow regime with zero av- 
erage flow in the a:-direction. At higher drives we find 
a crossover to a more ordered flow in Region IV, which 
can be identified by the reduction of fluctuations in Vy. 



2 



In region IV there can be a slight initial drift in the x 
direction due to the orientation of the moving particle 
lattice with respect to the triangular substrate, but as fd 
further increases, Vx goes to zero. 

In Fig. 2(c) we show Vx for the same system as in 
Fig. 2(a,b) but with slightly different initial particle po- 
sitions. Here the curve appears almost the same as in 
Fig. 2(b), but the flow in region II is in the positive x- 
direction. We find such symmetry breaking flow, which is 
equally likely to be in the positive or negative direction, 
for fillings N^/N > 1.0. There is a small region below 
the onset of region II where the flow shifts transiently 
between the positive and negative directions before it 
locks to one of the directions. We have done simulations 
for larger systems and find the same type of velocity- 
force curves; however, the transitory time for the flow 
to organize into region II increases with system size. At 
Nc/N = 2 and 3 and T = 0, region II does not occur be- 
cause the particle configurations are completely ordered 
at these commensurate fillings and there are no pertur- 
bations to knock the system into one of the symmetry 
directions of the lattice. At finite temperatures, however, 
region II reappears for these fillings. At incommensurate 
fillings there is a symmetry breaking due to the positional 
disorder of the particles. For fillings 1.0 < N/Nc < 2.0 
the velocity force curves look very similar to the curves 
shown in Fig. 2. For N/Nc > 2.0, region II can still 
occur; however, additional dynamical phases arise which 
appear as features in the velocity force curves and the 
noise fluctuations. We have also studied the tempera- 
ture dependence of the phases in Fig. 2. We find that 
there is a sharp transition temperature at which region 
II ends. This temperature is about O.TT^, where Td is the 
onset temperature for particle diffusion at fd — Q. In the 
inset of Fig. 4, we plot Vx for fixed Jd — 1.0 as a function 
of temperature showing the two branches of the velocity 
flow up to Tc where region II is lost. The sharp loss of 
region II as a function of T is consistent with the sharp 
transition from region II to the disordered region III. 

In Fig. 3 wc plot the particle positions and trajectories 
for the different phases. Fig. 3(a) shows the pinned phase 
for the system in Fig. 2, where a portion of the sites cap- 
ture two or more colloids. In Fig. 3(b) region II is shown 
for Fig. 2(b) at fd — 1-0, with the particles moving in ID 
paths along —30°. In Fig. 3(c), the case for Fig. 2(c) is 
shown for the same fd with the flow along +30°. Fig. 3(d) 
illustrates region III flow from Fig. 2(c) at fd = 3.2, 
where the particle positions and trajectories are disor- 
dered. Region IV flow at higher fd is much more ordered 
and the particles move mostly straight along the y di- 
rection; however, there are still some dislocations which 
cause the trajectories to be partially disordered. 

The region II flow can be regarded as an example of 
spontaneous symmetry breaking. There is no asymmetry 
across the .x-axis and in any given realization the system 
is equally likely to lock to either -|-30° or —30°. 




FIG. 3. Colloid positions (black dots) and trajectories 
(black lines) for the system in Fig. 2. (a) Pinned region I, 
/rf = 0. (b) -30° symmetry breaking flow at fd = 1- (c) +30° 
symmetry breaking flow at /d = 1. (d) Disordered flow in 
region III, fd = 3.2. 



Spontaneous symmetry breaking in general is not limited 
to equilibrium systems. In condensed matter, the most 
familiar example is that of a ferromagnet in zero external 
field where there are two energy minima. A small pertur- 
bation from temperature or dynamics biases the system 
in one direction or another. In our system, when the flow 
follows one of the symmetry directions, the flow is more 
ordered and the particle spacing is more uniform. 

We have also considered a single particle for the system 
in Fig. 1 as well as other fillings with Nc/N < 1.0. In 
these cases the global flow locking does not occur. This 
is due to the lack of interactions at these low fillings. At 
depinning, each particle can move in cither the positive 
or negative direction at each potential maximum. If there 
are no neighboring particles blocking one of the routes, 
the particle flows in a random zig-zag pattern. 

Some of the features in the velocity force curves in 
Fig. 2 can be understood by force balance arguments. 
For the commensurate case Nc = N, the colloid config- 
urations are completely symmetric and the interaction 
forces are zero, so depinning occurs at fd = /™°^) where 
jmax _ g g jg ^^ic maximum force force from the sub- 
strate. For higher Nc, some substrate minima capture 
more than one colloid. In these dimerised phases, a col- 
loid feels an additional force from the other colloid in the 
same minima. This will reduce the force needed to depin 
the colloids by the interaction force fij{ao/2). Using our 
parameters this would give a depinning force of 0.8 which 
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Time 

FIG. 4. The transverse velocity 14 for the ramp up phase 
of 0.0 < fd < 4.8 for time < t/U < 1; down from fd = 4.8 
to 0, 1 < t/tr < 2; and up again to fd = 4.8, 2 < t/tr < 3. 
Here tr = 2.24 x lO'' molecular dynamics time steps. Inset: 
transverse velocity for fixed fd = 1.0 vs T/Td. 

is somewhat higher than the measured value due to the 
collective interactions of many dimers. Fig. 3(a) shows 
that the dimcriscd states preferably orient cither in the 
positive or negative 30° direction, so that is the direction 
in which they depin. The transition to region IV flow 
occurs for fd > /J""^ where the colloids can move freely 
over the substrate potential maxima. 

We next consider the hysteretic properties of the 
phases by performing simulations where the external 
drive is cycled up, down, and up again. In Fig. 4 we 
show the results for the same system as in Fig. 2 for 
three drive cycles. For < t/tr < 1 where t,. = 2.24 x 10^ 
molecular dynamics time steps, fd is increased from zero 
to fd = 4.8. On the ramp down, beginning at t/tr = 1) 
there is clear hysteresis in the III-II transition with region 
III persisting longer on the downward sweep. Hystere- 
sis across the II-III transition is consistent with the first 
order-like sharp jump in Vx across the II-III transition. 
On the ramp down, region II is again in the negative 
direction. The slope is also smaller as there are more 
commensurate pinned rows. We note that because the 
system goes through region III the flow is disordered so 
that the pinned phase near the end of the ramp down 
phase is not necessarily the same as the initial pinned 
phase. On the next ramp up, beginning at t/tr = 2, Vx 
in region II reverses and is in the positive direction. For 
continued cycling, even at T = 0, the direction of the 
region II flow falls randomly in either direction. This is 
due to the dynamical disordering introduced by the re- 
gion III flow, which destroys the memory of the previous 
pinned phase. We observe similar hysteresis effects for 
other flUings in the range 1 < Nc/N < 4. 

In addition to the colloidal system, essentially the same 



situation occurs for vortices driven in superconductors 
with periodic pinning arrays. In this case we considered 
each pinning site to capture one vortex. At fields with 
more than two vortices per pinning site, the same dimer 
states occur for vortices in the interstitial regions between 
the pinned vortices and the same phases arise. 

To summarize, we have found a novel spontaneous 
symmetry breaking phenomena for collectively interact- 
ing driven particles on triangular substrates. The global 
particle flow occiirs along one of two symmetry directions 
and not in the direction of the applied drive. This ap- 
pears as a positive and a negative finite Hall velocity. The 
flow along the symmetry directions is much more ordered 
than the flow in the direction of the drive. We also find 
a series of dynamical phase transitions between ordered 
and disordered flows which produce jumps in the longi- 
tudinal and transverse velocity force curves. The sym- 
metry breaking combined with the flow transverse to the 
drive produces interesting hysteresis phenomena in the 
transverse velocity curves where the particles may lock 
to different directions upon increasing and decreasing the 
drive, as long as the particles pass through a disordered 
flow regime. The results in this work can be realized for 
colloids driven over optical trap arrays and vortices in 
superconductors with periodic pinning arrays. Similar 
effects may be possible for atomic friction systems. 
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